Multiple glass transitions in star polymer mixtures: 
Insights from theory and simulations 
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The glass transition in binary mixtures of star polymers is studied by mode coupling theory and 
extensive molecular dynamics computer simulations. In particular, we have explored vitrification in 
the parameter space of size asymmetry 5 and concentration pi of the small star polymers at fixed 
concentration of the large ones. Depending on the choice of parameters, three different glassy states 
are identified: a single glass of big polymers at low S and low p2, a double glass at high S and low 
p2, and a novel double glass at high p2 and high 5 which is characterized by a strong localization 
of the small particles. At low S and high p2 there is a competition between vitrification and phase 
separation. Centered in the (5, p2)-plane, a liquid lake shows up revealing reentrant glass formation. 
We compare the behavior of the dynamical density correlators with the predictions of the theory 
and find remarkable agreement between the two. 



I. INTRODUCTION 

The slowing down of the dynamics in a supercooled liq- 
uid near the kinetic glass transition still requires a com- 
plete molecular understanding though recent progress 
has been made by experiments [U [21 [3] , by computer sim- 
ulations [1 M EE El [H H] and theory [M EEU H EH ■ The 
inherent polydispersity of colloidal systems facilitates the 
formation of disordered arrested states, as compared to 
the formation of the underlying equilibrium crystal [lj. In 
simulations, it is difficult to vitrify a pure one-component 
system of spherical particles due to the onset of crys- 
tallization, so that binary mixtures constitute a much 
richer and more convenient system to study the glass 
transition [T3] . The thermodynamically stable state of an 
homogeneous undercooled binary melt typically involves 
two phase-separated binary crystals, A-rich and B-rich 
[151 116] or a complex alloy in which the two components 
are mixed [23 [T7] . In either case, the processes involved 
to reach true thermodynamic equilibrium require collec- 
tive, mutually coordinated motions that take enormously 
long time. Therefore, the glass transition does take place 
also in simulation, since crystallization is prevented by 
the inherent polydispersity of the mixture. Hence, either 
binary or polydisperse systems have always been used as 
prototype models for the glass transition, for example re- 
pulsive soft spheres [U [TSJ QI5] or Lennard- Jones mixtures 
[201 EH E2], charged mixtures [23], hard spheres [24], and 
Yukawa systems [251 US] ■ I n fact, the main reason for con- 
sidering these slightly asymmetric binary mixtures is that 
crystallization was avoided. Upon cooling down typically 
all species freeze into a glass at the same temperature. 

More recently, insight has been gained into the glass 
transition of strongly asymmetric mixtures consisting of 



large colloidal spheres and small additives such as smaller 
colloids or non-adsorbing polymers which induce an ef- 
fective depletion attraction between the large particles 
whose range is dictated by the size of the small parti- 
cles. For strong asymmetries, very short range attrac- 
tions (relative to the core size of the large particles) can 
thus be realized. A new type of glass was observed for 
these attractive interactions — termed an attractive glass 
— in experiments of colloid-polymer mixtures [271 and in 
computer simulations [25] , following predictions based on 
mode coupling theory [29] , There is an interesting reen- 
trant behavior [30 for increasing concentration of the 
additives: the repulsive glass (at vanishing additive con- 
centration, i.e. in the absence of any attraction) melts 
upon a threshold additive concentration. The resulting 
melt re vitrifies into a second type of glass, namely the 
attractive glass, for higher additive concentrations. Both 
for repulsive and attractive glass, only the large particles 
freeze into a disordered glassy matrix while the small ad- 
ditives always remain mobile and stay in a fluid state oc- 
cupying the "holes" within the glassy matrix. This has 
been made explicit by computer simulations and mode 
coupling theory for the Asakura-Oosawa model in which 
the polymer-polymer interaction is assumed to be ideal. 
It was found [3T] that the bare mobility of the additives 
is very important: when it is comparable to that of the 
colloids there is a freezing into a common glassy matrix 
of colloids and additives, while when it is much larger 
than the one of the colloids, freezing involves only the 
colloid component. 

Experiments on binary hard sphere mixtures 32J have 
shown that there is a geometrically-inspired cross-over 
behavior in the dynamics of the small spheres when the 
ratio 5 of the two spheres diameter is about 0.15. This is 
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well below the size asymmetry considered in earlier mode- 
coupling calculations [33] - Another system in which dy- 
namical arrest has been studied theoretically and ex- 
perimentally are star polymers and micellar solutions 
[331 US US H3- The first systematic attempt to vary 
the size asymmetry and to monitor its influence on the 
glass transition was performed in a binary mixture of star 
polymers in a good solvent [38 . Mode coupling theory 
and experimental data were compared for different size 
ratios 5 of the star polymer species. At low S the small 
additives remain ergodic in the glassy matrix formed by 
the large stars (single glass) while at high 5 both species 
become arrested simultaneously (double glass). 

For increasing concentration of the smaller star poly- 
mers a melting into a fluid was found. The threshold 
concentration of the smaller stars for melting behaves 
non-monotonically as a function of S for fixed concen- 
tration of the big stars with a [/-type shape. This non- 
monotonicity is a result of the ultrasoft effective interac- 
tions between the star polymers [39] and has no coun- 
terpart in hard-sphere based systems. Additive induced 
melting has also been found in mixtures of star poly- 
mers and liner polymer chains [401 141] . More recently, a 
detailed study of the structural properties at high densi- 
ties of the additives was performed [42] and an additional 
double glass, called asymmetric glass, has been identified. 

In the present paper we elaborate more on the sys- 
tem of binary star polymers. Our motivation to do so is 
twofold: Firstly, aiming at a deeper understanding of the 
full dynamical behavior of collective- and self-correlations 
close to the various glass transitions, we performed ex- 
tensive molecular dynamics computer simulations and 
compared the results with predictions of Mode Coupling 
Theory calculations. In this fashion, we were able to 
extract the dynamical scaling laws by fitting the simula- 
tion data. As regards the arrested plateau values of the 
dynamical correlations and localization lengths of both 
large and small components we find quantitative agree- 
ment between simulations and MCT. Secondly, we have 
expanded the parameter space up to concentrations of 
added small stars higher than previously studied [35J and 
found a reentrant vitrification of the molten state [32 . In 
the plane of size asymmetry S and added concentration 
P2 the molten state formes a lake with an O-shape. In 
addition to the asymmetric glass, we also observe an at- 
tractive glass at low S and high pi , which is unstable with 
respect to phase separation. 

The paper is organized as follows: in Sec. II we intro- 
duce the theoretical methods used in this paper. Sec. 
Ill contains the results based on Mode Coupling Theory, 
where we distinguish four regions in the (5, p 2 ) plane. 
Simulation data for the glass transition line based on iso- 
diffusivity curves are presented in Sec. IV, whereas sim- 
ulation data for the dynamical correlations in the four 
different regions arc subsequently discussed in Sec. V. 
Finally, in Sec. VI we summarize and present our conclu- 
sions. 



II. THEORETICAL METHODS 

We describe the system in terms of effective interac- 
tions between the star centers [39] . The monomer degrees 
of freedom have been traced out in this approach. This 
method is applicable to the problem at hand, because we 
are interested in the arrest of the stars as whole objects 
and not in the arrest of the monomers themeslves, which 
occurs at much higher concentrations. Concomitantly, 
the monomers maintain their fluctuating nature through- 
out. However, possible entanglements of the arms are not 
included in this description. The effective interaction be- 
tween star centers depends parametrically on the num- 
ber of arms (or functionality) /, and shows a logarith- 
mic divergence when the two stars overlap within their 
corona diameters cr-j, i — 1,2, (where a = Rh [40] , with 
Rh the hydrodynamic radius), followed by a Yukawa-tail 
for larger separations [39"] . 

In this work we study binary star mixtures with func- 
tionalities fx — 263 and f 2 = 64. We express the mixture 
composition by quoting the values p\a\ of the big and 



P2&1 of the small stars, where pi 



1, 2 are the respec- 



tive number densities. The size ratio S = <y%l<J\ is an ad- 
ditional parameter of the mixture. The two-component 
generalization of the aforementioned effective star-star 
potential reads [43] . 
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In order to calculate the fluid structure of the mixture 
we solve the binary Ornstein-Zernike (OZ) equation [44] 
with the Rogers- Young (RY) closure [35]. This amounts 
to postulating the additional relation to hold between the 
total correlation functions hij(r) and their direct coun- 
terparts Cij(r): 



h ij (r) = exp[-(3V ij (r)} 
exp[(%(r) 



-^(r))/(r)]-l 
f(r) 



- 1, 



(3) 



where f(r) — 1 — exp(— ar) with (a > 0). The value of 
the parameter a is determined by requiring the equality 
of the 'fluctuation' and 'virial' total compressibilities of 
the system [35] [3S]. For a — > the Percus-Yevick closure 
is obtained, while for a — > oo the hypernetted chain clo- 
sure (HNC) is recovered [44]. The RY closure has been 
shown to correctly describe liquid properties for a num- 
ber of repulsive potentials, including star polymers [47J, 
spherical ramp [48 and square shoulder [49] potentials. 
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From solving the OZ equation, it is also possible to 
calculate the partial static structure factors, describing 
the density correlators at equal time, 



Sij(q) = S tj + y/PiPjhij(q) = 



1 



(p*(q,0) ft (q,0)) 
(4) 



where hij(q) denotes the Fourier transform of the total 
correlation function and 



Pi 



(q,i) =^exp iq-r, W) (t) 



1=1 



with rj being the coordinates of the Z-th particle of 
species j (j = 1,2) and the asterisk denoting the complex 
conjugate. The partial static structure factors calculated 
within RY approximation for binary star polymer mix- 
tures have been found in very good agreement with those 
obtained by simulations employing the effective potential 
of Eq. ffl) [50]. 

Our theoretical study of the glass transition is based 
on the ideal Mode Coupling Theory (MCT), a theory 
that describes the time evolution of the density auto- 
correlation functions, starting only from the knowledge 
of the static structure factors, via a set of coupled 
integro-differential equation [TO] . The normalized time- 
dependent collective density autocorrelation functions 
are defined as, 

^(g,t) = (p|(q ) 0) ft (q,t))/(p|(q ) 0) ft (q,0)). (5) 

It is also useful to focus on the self part of the density 
correlation functions, which describes the dynamics of a 
tagged particle, 
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The long-time limit values of (f>ij(q,t), i.e. the partial 
non-ergodicity factors, respectively for the collective and 
self correlation functions, are defined as, 



t — t-OO 

f?(q) = lim ^(q,t). 



(7) 
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A glass transition is identified within MCT as an ergodic 
to non-ergodic transition, when the non-ergodicity factor 
discontinuously jumps from zero, typical of a fluid, to a 
finite value typical of a glass [TU] . 

In binary mixtures with large size asymmetries, there is 
a vast separation of time-scales between large (species 1) 
and small (species 2) stars due to the enormous difference 
in mobilities. Therefore we can use one-component MCT 
for the large stars only, and solve the equation for the 
non-ergodicity parameter f\\(q) |31j . 
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V(q,k)/u(fc)Ai(|q-k|), (9) 



with 



V(q,k) = ^[q-(q-k)cii(|q-k|)+q.kcn(A;)] : ' 



x 5n(g)5n(fc)5ii(|q-k|), 



(10) 



The solution with the largest fn(q) is the long-time limit 
[TU] of the density correlators. We note here that, in the 
one-component effective treatment, we use Sn(q) calcu- 
lated from solving the binary OZ equation within RY 
closure, arising from the full binary interactions of Eq. 
([lj. In this respect, we are treating the small stars as an 
effective medium, and their influence on the interactions 
between the large stars is taken into account explicitly 
and not through an effective one-component picture [46] . 

MCT assumes a vast separation of time scales between 
the slow and fast degrees of freedom. For small size 
asymmetries, the separation of timescales between the 
different star species decreases and one-component MCT 
cannot be used |31j . In the general case of a mixture, 
when relaxation time-scales of the two species are com- 
parable, the long-time MCT equations (Eq. |9| are easily 
generalized [19j [5TJ in terms of matrices as, 



f(q)=S(q)-{S(q)- 1 +H[f],q)Y 
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where S(q) — ^/XiXjSij(q), with Xi = ^/(TVi+A^) being 
the number concentration of species i. 
The functional J-( [f ] , q) is defined as 
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with p = P1+P2 andp = |q— k|. The vertices I^ mm '(q,k) 
depend on the equilibrium structure of the system and 
are given by 
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In the case of binary or multicomponent MCT, the the- 
ory predicts, through the couplings between species in 



Eq. (Ilia simultaneous jump from zero to a finite value 
for all partial collective non-ergodicity parameters fij(q)- 
However, this does not hold for the self part fi(q), for 
which the equations are |52j . 
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so that it is possible to distinguish also the case of mobile 
particles in a frozen environment, for example for the 
small stars in our case, when f2i(q) 7^ but = 



In general multiple glassy solutions of Eq. ( 12 1 can ex- 



ist, due to higher order bifurcation singularities. This 
case has been extensively studied for effective one- 
component systems of attractive colloids, where higher 
order MCT singularities have been predicted theoreti- 
cally |54j and verified by numerical simulations |55[ 156] 
and experiments [S7J US [35] . For multicomponent sys- 
tems, calculations of MCT higher order singularities have 
not been reported so far. However, hints of the presence 
of such higher order singularities and glass-glass transi- 
tions, in particular of one of the involved species, have 
emerged in simulations of asymmetric polymer blends 
and soft sphere mixtures [59] [60] [61]. In the case of 
distinct glassy states, the non-ergodic properties of the 
various arrested states will be different from one another, 
a feature that will be reflected on the corresponding non- 
ergodicity factors fij(q,t). Also, close to such singulari- 
ties, competition between the different glassy states mod- 
ifies the standard MCT dynamical behavior close to the 
liquid-glass transition (labeled usually of type B or A 2 ). 
For such an A 2 transition, a two-step relaxation is well 
described by MCT, through an asymptotic study of the 
correlators near the ideal glass solutions. 

For a generic density correlator 4> in the liquid the crit- 
ical non-ergodicity parameters f c (q) can be extracted. 
The departure from the plateau, i.e. the start of the 
a-process, is described by a power law, called the von 
Schweidler law, regulated by the characteristic exponent 
b, which is independent of the particular q- vector consid- 
ered, 

4>(q,t) - f C (q) ~ -h^(q)(t/r ) b + (q) (t / r ) 2b (15) 

with To being the characteristic time of the relaxation. 
The quantities h^(q) and h^ 2 \q) are referred to respec- 
tively as critical amplitude and correction amplitude 1 OJ . 
The a— relaxation process can also be well described by 
a stretched exponential, i.e. 



<j>(q,t) = A(q) exp[-(i/T,y 



(16) 



where the amplitude A(q) determines the plateau value, 
and the stretching exponent j3 q < 1. For large q values, 



it has been shown that (3 q 



The von Schwei- 



dler exponent b is related to another important MCT 
parameter, named the exponent parameter A 10 . This 
includes all information on the dynamics of the systems 
on approaching the MCT transition. It takes values in 
the interval 1/2 < A < 1. 

Close to higher order singularities, the dynamical be- 
havior is predicted to obey a logarithmic behavior and A 
should approach 1. In particular, it holds, 



4>{q, t) ~ f c (q) ~ Kq) Lb« ln(i/r) + B< 2 > ln 2 (i/r) 



We define q* as the wave-number at which the correla- 

(2) 

tors decay purely logarithmically, that is B q , = 0. Here 
r stands for a time-scale which diverges if the state ap- 
proaches the mathematical singularity. The formula is 
obtained by asymptotic solution of the MCT equations 
[63] . The first term f c (q) is the sum of the non-ergodicity 
parameter at the singularity plus a correction which de- 
pends on the distance from the singularity. 

We will use Eq. (15l-([l7| to describe the the corre- 
lators (f>n(q,t) and (f> 22 (q, t) independently. In the fol- 
lowing, for those states where the density correlators do 
not follow power-laws or stretched exponentials, we adopt 



the expression in Eq. (17) to extract the non-ergodicity 



(17) 



parameters from the simulations and find that they are 
indeed very well described by that function. 



III. MCT RESULTS 

Following the experimental setup discussed in |38j . we 
investigate the effect of the addition of a second com- 
ponent of varying size and concentration on the glass 
formed by the one-component large stars with /1 = 263. 
As we have shown before [35], the qualitative trends do 
not depend on details of the mixture properties, there- 
fore we consider only f 2 — 64 in this work. We com- 
mence from a glass of large stars only, with a fixed den- 
sity fi\o\ = 0.345, corresponding to a glassy state within 
MCT [34 . At this state point, we add small stars and 
study the stability of the large star glass [75 . 

From the MCT analysis, we determine a glass transi- 
tion diagram in the (S, p 2 ) plane (see Figure [lj. Results 
obtained both in the full binary treatment of the mix- 
ture and in the effective one-component representation 
are shown. With respect to the (p, T) phase diagrams 
that are usually shown in colloid-polymer (CP) mixtures 
[27] . we focus here on a single isochore and consider the 
addition of small stars of different sizes and concentra- 
tions. 

Examining Fig. [l] a liquid region is found, being sur- 
rounded by many different glasses. We detect a reentrant 
behavior on increasing the density of additives both for 
small and large sizes of additives as well as a reentrance 
in size ratio. For convenience, we can roughly divide the 
phase diagram of star-star mixtures in four different re- 
gions to be discussed in the following: 

1 . small size and small concentration of the additives 
(region I): single glass; 

2. large size and small concentration of the additives 
(region II): double glass; 

3. large size and large concentration of the additives 
(region III): asymmetric glass; 

4. small size and large concentration of the additives 
(region IV): attractive glass/phase separation. 
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FIG. 1: Ideal MCT kinetic phase diagram in the size ratio- 
additive density plane. Symbols are the bracketing state 
points for the ideal-glass transition, determined by using the 
two-component MCT. The full line through the points is a 
guide to the eye. The dashed line is the ideal-glass line calcu- 
lated on the basis of the one-component MCT. Four different 
regions with glassy states can be identified: we name them 
respectively single glass (region I), double glass (region II), 
asymmetric glass (region III) and attractive glass/phase sep- 
arated glass (region IV). These calculations refer to the case 
piaf = 0.345. The dotted lines roughly separate the different 
regions of the phase diagram. 



Preliminary information on the single, double and 
asymmetric glasses have been reported in Ref. [3EJ HI] ■ 
Here we summarize the most important aspects for com- 
pleteness. We start by examining region I. The small 
stars always remain liquid in this regime, even in those 
states in which the large ones are arrested (single glass) 
|38j . The addition of small stars with size ratio up to 
S w 0.25 causes the large glass star to melt due to a 
depletion-induced softening of the repulsions with respect 
to the one-component case [JBJ, weakening the stability 
of the cages [35] . 

On the other hand, moving into region II by increasing 
the small stars diameter, the addition of a second compo- 
nent of stars with 0.25 < 5 < 0.45 leads to a glassy state 
formed by both components (double glass) [35]. Here, 
an increase of p2 also melts the glass but through an en- 
tirely different mechanism, found in MCT and also con- 
firmed by experiments [38 . Indeed, the glass becomes 
stiffcr through the addition of a second glassy compo- 
nent, as shown by the growth of the clastic modulus in 
theory and experiments reported in |381 150] up to the 
melting density. At this point, it becomes entropically 
more convenient to free available volume for both stars 
simultaneously, which amounts to a finite energetic in- 
crease, sufficient to cause liquification. Therefore, this 
melting mechanism is only possible through the softness 
of interactions involved, and it does not have an equiva- 
lent in terms of hard-sphere mixtures. We note, that the 
presence of two different glasses for asymmetric mixtures, 
a single to a double glass, depending on the size ratio, was 



already found for hard spheres |55J |55] , but the new 
feature of the presently investigated soft mixtures is that 
both glasses can be simultaneously melted by increasing 
the amount of the small additives in the mixture. 

Clearly, a theoretical framework capable to correctly 
describe the single glass should be based on the slowing 
down of the large stars only, i.e. one-component MCT, 
where the smaller stars are assumed to form a fluid 
medium that does not participate in the glass formation, 
is the conceptually right framework to use. On the other 
hand, in the double glass both species become compa- 
rably slow, thus compelling the use of binary MCT. It 
turns out that in the case of star polymers mixtures, dif- 
ferently from that of colloid-polymer mixtures |31j . at 
low 5, i.e., in the single glass regime, both treatments 
give qualitatively the same behavior, as illustrated also 
in Figure [T] while at large 5, in the double glass regime, 
the discrepancy becomes large, signaling the tendency of 
the small component to also arrest. Only binary MCT 
can be used in this region, giving rise to the U-shape also 
observed in the experiments for low p 2 [38j . Hence, in 
the following we just refer to two-component MCT in all 
regions, unless explicitly mentioned. 

In the present work we complement previous results 
[3"51 150] with the investigation of Regions III and IV. 
In region III, for large concentration of additives with 
6 ~ 0.4, we find a reentrant vitrification with respect 
to p2, with the appearance of a new glassy state with 
smaller localization length for large 6 and large p2- Due 
to the softness of the interactions local rearrangements 
take place which leads to a change of coordination of the 
large stars as well as to an asymmetric deformation of the 
cages )42j . The increase of p2 causes a second jamming 
at a new length scale that will depend on the particular 
parameters of the mixture. 

Finally in region IV, at low S, the depletion mechanism 
induces effective attractions between the large stars [46] . 
causing phase separation. Phase separation is accom- 
panied by an increase of the structure factor at infinite 
wavelength, and a blind use of it in the mode-coupling 
equations would predict glass formation driven by the 
small wavevectors. However, with the use of RY struc- 
ture factors, we are only able to find a glass within binary 
MCT for S = 0.1, while the one-component treatment 
always provides a liquid state up to the maximum inves- 
tigated P2 approaching phase separation. We could not 
extend the analysis to larger p2 due to a breakdown in 
the convergence of the Rogers- Young closure, as well as 
the simplest HNC closure in this region, a possible sign 
of a close thermodynamic instability. 

In addition to the phase diagram, we report in Fig. [2] 
the calculated MCT partial non-ergodicity parameters, 
both for self [Fig. [2^a-b)] and collective [Fig. [2jc-d)] re- 
laxation dynamics for 5 = 0.1 and S — 0.4. Such results 
correspond to the glassy points closest to the transition, 
hence they are 'critical' non-ergodicity parameters. The 
small stars /|(g) allows us to distinguish immediately 
a single from a double glass, since it is identically zero 
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FIG. 2: Partial non-ergodicity parameters at the MCT tran- 
sition. The (a) and (b) panels show respectively the large-star 
and small star components of the self non-ergodicity param- 
eters for the three different glasses. They also show Gaussian 
fits (dotted lines), which allow us to extract the localization 
length. The f( of the large stars in the single and double 
glass are virtually identical, therefore the lines fall on top of 
each other, while those for the small stars are very different 
in all three cases. In the (c) and (d) panels we show the 
non-ergodicity parameters for the collective dynamics. For 
all curves pio-f = 0.345. For the single glass P20"i = 0.12 and 
8 — 0.1, while in the double glass p2<J\ = 0.05 and 8 = 0.4. 
The asymmetric glass corresponds to p2<jf = 1.27 and 8 — 0.4. 



for the asymmetric glass, i.e., the cage size reduces by 
a factor of two. A larger variation in Zmct is predicted 
for the small stars from 0.51cti for the double glass at 
low pi to O.llcri for the asymmetric glass at high p2, a 
value which is almost identical to the one for the large 
stars in the same state. For the collective fu(q), the re- 
sults are very similar to those of the self non-ergodicity 
parameters, only much stronger oscillation peaks are ob- 
served in phase with the static structure factor. For the 
small stars, these appear only in the asymmetric glass 
regime. We notice that upon changes in S within the low 
P2 glasses, we detect no change in fu(q) between single 
and double glass. However, a significant increase is ob- 
served for large 5, upon increasing p2, while approaching 
the asymmetric glass. This increase results in a more 
extended q- width in /n(q), a finding similar to that for 
attractive glasses in colloid-polymer mixtures, where fur- 
ther increase of the additives leads to an increase of the 
non-ergodicity parameter both in amplitude and range 
|28j . Therefore, in the asymmetric regime, a strong lo- 
calizing mechanism is at hand, with some similarities to 
an effective attraction. However, looking at the partial 
Sij(q) we do not see an increase in the compressibilities 
or a non-monotonic dependence in the growth of peaks. 
Simply, the structure factor peaks for large stars continue 
to decrease upon addition of small stars, while these dis- 
play themselves more and more enhanced peaks. In |42J, 
we have shown that the change of structure is accompa- 
nied by a dramatic decrease in the coordination number 
of the large stars and an asymmetric deformation of the 
large star cages. 

Finally, an analysis of the MCT exponents reveals that, 
for both single and double glass, A ~ 0.71 and b ~ 0.625, 
the same result obtained for one-component star polymer 
solutions. However, we note that A grows considerably 
close to the asymmetric glass transition. At the transi- 
tion point corresponding to 6 = 0.4, A ~ 0.8, correspond- 
ing to a von Schweidler exponent b ~ 0.47. 



IV. SIMULATIONS: NUMERICAL DETAILS 



for the single glass, as already found in binary asym- 
metric hard sphere mixtures |52j . Thus, although, in 
binary MCT the broken ergodicity transition is postu- 
lated to happen simultaneously for both species, the self 
non-ergodicity parameter allows to distinguish between 
mobile and immobile species. Also, the partial f£(q) for 
both types of stars have a larger width in the asymmetric 
glass than in the double glass. Fitting f?(q) to a Gaus- 
sian curve exp(— <7 /6ZjJjct) P rov ides an estimate for the 
localization length Imgt within the glass. Fits are also 
reported in the figure. The Gaussian fit works very well 
for large stars, but deviations are observed for small stars, 
in both double and aymmetric glasses. The extracted lo- 
calization length is Imgt ~ 0.26<7i for large stars both 
in single and double glass, while it is Imgt ~ 0.1 4ci 



We performed standard Newtonian dynamics simula- 
tions for stars interacting via the effective potential in 
Eq. ([IJ . Although the use of Brownian dynamics simu- 
lations would be more realistic for star-star mixtures to 
mimic the solvent effects, this is less efficient and it would 
not allow us to probe the full slow dynamics regime effi- 
ciently. Moreover, several confirmations of independence 
of long-time properties from microscopic dynamics have 
been provided so far [SH ES] for systems close to their 
glass line. We investigate several compositions, in order 
to have a complete picture of the phase diagram and of 
the dynamical slowing down. Due to the various compo- 
sitions to be investigated, the number of particles varies 
throughout the diagram. In general, we fix N\ = 1000 
and vary N 2 accordingly. However, for pi<j\ > 10, we use 
Ni = 250. Also, for p 2 (jf < 0.1, i.e. in single and dou- 
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ble glass regimes, we fix the minimum number of small 
particles to N 2 = 500 and vary N\ accordingly. We in- 
vestigate fi — 263 and f 2 — 64, the same values used 
in the MCT analysis. To improve statistics, we average 
over 5 independent runs for each studied state point. 

In order to describe the dynamics of the system prop- 
erly, we have to choose the short-time mobilities of the 
small particles with respect to the large ones. Follow- 
ing the ideas in [3T], we use the scaling for the mass 
to typical of the star polymers, namely to ~ J 2 / 3 ^ 5 / 3 
[SB"] 167] , This results in a ratio between the small and 
large particles mass equal to m 2 /m\ = (/2//i) 2 ^ 3 ^ 5 ^ 3 - 
Note that a physical choice of the small particle mass 
is crucial, since as already noticed in earlier simulations 
of explicit colloid-polymer mixtures [31] . the addition 
of small particles of mass identical to the one of the 
large particles, has the result of further slowing down 
the large particle dynamics. This would imply no melt- 
ing at any given p 2 and 5 in contrast with MCT predic- 
tions and experimental results. Units of length, mass 
and energy are <7i, mi and k^T, respectively. Time 
is measured in tmd = / (k^T) and the integra- 

tion time step is varied according to the mass ratio, as 
5t = 5 • 10 _3 v /to 2 /to 1 t md . 

We employ numerical simulations to (i) evaluate the 
iso-diffusivity lines in the p 2 — S plane and (ii) study 
specific features of the dynamics for selected state points 
approaching the four glass regions. 

Simulations for the isodiffusivity lines are performed 
for p\o\ = 0.345. Note however that at this density the 
slowing down of the dynamics observed in the simula- 
tions is not very significant. To study the dynamics on 
approaching the glass transition(s) we have been forced 
to increase the density of the large particles to generate 
a slower dynamics. Indeed, similarly to other studied 
cases like that of hard spheres, a systematic shift of the 
MCT transition with respect to the slowing down taking 
place in experiments or simulations, is present. This shift 
usually applies in the direction of MCT anticipating the 
onset of glassy dynamics and our findings agree with this 
general picture. 

Moreover, we have considered a small poly-dispersity 
in the large particles diameter. This is particularly im- 
portant at low p 2l when the amount of small stars is 
not sufficient by itself to prevent crystallization of the 
large stars. We selected the value p\o\ = 0.41 to in- 
vestigate the dynamics of single and double glass respec- 
tively, while we retained the value p\o\ = 0.345 for the 
asymmetric glass. The larger density and the absence of 
crystallization allows us to observe a slowing down of the 
dynamics for the large particles (even in the absence of 
additives), and consequently also the single and double 
glass. We employ a polydispersity in the large-stars di- 
ameter, chosen to follow a Gaussian distribution with a 
width of 10% [M]- 

We have checked that, compared to the one-component 
case, the dynamics are not affected. For the compar- 
ison, we have to account for the polydispersity in the 
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FIG. 3: Diffusivity of the one-component system as a func- 
tion of star density. Circles correspond to the monodisperse 
system with / = 263 and squares to simulations with added 
polydispersity to prevent the system from forming a crystal. 
The exploration of three orders of magnitude in D\ is thus 
possible. 



right way when calculating the density, i.e. by matching 
the value p\ (of) with the corresponding one-component 
density. In this way, we can decrease the diffusivity by 
three orders of magnitude by increasing p\ , as illustrated 
in Figure [3j comparing results for the monodisperse and 
polydisperse one-component reference system. 

Finally, for the simulations of the attractive glass, we 
have focused on monodisperse large stars, due to the ab- 
sence of systematic crystallization in this case. 



V. SIMULATION: ISODIFFUSIVITY LINES 

The iso-diffusivity lines are precursors of the glass 
transition lines, and they always display a similarity in 
shape as the MCT line, as demonstrated already for one- 
component star polymers [34], for other repulsive poten- 
tials |35], as well as for short-range attractive square well 
potentials [68 . 

To study iso-diffusivity lines we calculate the mean- 
squared displacement (MSD) (rf(t)) (i = 1,2) of large 
and small stars separately. From the long-time be- 
havior, we extract the respective self-diffusion coeffi- 
cient defined as, following the Einstein relation, Di = 
lim t ^ o(r?(t)>/(6t). 

In Fig. [4] we report the iso-diffusivity lines for three 
different values of Di/Dq in the (6, p 2 ) plane, where Do 
is the diffusion coefficient in the absence of small parti- 
cles. The fact that we examine only a decade in decrease 
of Di/Dq is due to the many constraints to study the 
whole phase diagram under very different mixture com- 
positions. The line corresponding to the slowest states 
is drawn only on the large 5 side, because the diffusivity 
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FIG. 4: Iso-diffusivity lines from MD simulations in the 
(p2, ^-representation. The density of large stars is kept con- 
stant at piaf — 0.345 and the large stars are monodisperse. 
The line for the smallest value of the diffusivity is not closed 
because the chosen value of D\ is smaller than the diffusiv- 
ity of the reference one-component system (Do)- The dashed 
curve is an isodiffusivity line of the small stars. 




FIG. 5: Mean squared displacements for large and small stars 
in the single glass. The densities are piaf — 0.41 and p2o\ = 
0.1; the size ratio is 5 — 0.1. The horizontal line indicates 
the estimate of the squared localization length for large stars, 
defined as the inflection point of the MSD. 



is lower than in the one-component system and therefore 
no closed loop exists for this value of D\. Notwithstand- 
ing these details, the emergence of an asymmetric "O" , 
that resembles the shape of the MCT phase diagram, 
is evident. This was confirmed also for higher p\ [42], 
demonstrating that the topology does not depend on the 
precise value of p\ studied. Recently, this asymmetric 
"O" has been observed experimentally |4"2"] . 

Additionally, from the simulations we can also trace 
out correctly the iso-diffusivity lines of the small stars 
(Fig. [4|. Within binary MCT, there should be no dis- 
tinction between the glass lines of the large and small 
stars, based on the assumption that small stars are as 
slow as the large stars. It has a similar shape as the 
large <5/large pi region. Note that the dynamics of the 
small stars is much faster for low 5, where the mobility 
of the two species is very dissimilar, and hence the small 
star isodiffusivity lines do not show an "O" shape. 



VI. SIMULATION: MULTIPLE GLASSY STATES 

We can now exploit the simulation results to directly 
probe the four regions of slowing down and compare the 
different glasses, in order to clarify the characteristics 
of the multiple glassy states appearing in these types of 
mixtures, in comparison with the MCT results. 



A. Region I: Single glass 

We consider the state point pi = 0.41, p 2 = 0.1 and 
5 = 0.1 (2250 particles in the simulation box). The 



mass ratio is TO2/7711 = 8.4 x 10~ 3 , allowing an effec- 
tive one-component description [31 . Indeed, the mean 
squared displacement in Fig. [5] shows a clear separation 
in time-scales between the two species. After the ini- 
tial ballistic regime, the small stars simply become diffu- 
sive, while the large ones display a clear plateau, lasting 
about three decades in time, before eventually reaching 
the diffusive regime. Correspondingly, a difference in the 
diffusion coefficient (and relaxation time) of about three 
orders of magnitude is observed. By convention, we de- 
fine the localization length as Imd — \J { r2 (t*)), where t* 
is the point of inflection of the MSD in the log — log plot. 
We therefore provide evidence that the large stars are 
nearly arrested in a jammed state, with a localization 
length that can be extracted from the plateau height 
!md ~ 0.25(7i (horizontal line in Figure [5]). A similar 
value could be extracted from a simulation of star poly- 
mers in theta-solvent [64]. Also our MCT results dis- 
cussed provide a similar value, suggesting rather robustly 
that the localization length of a soft star glass is slightly 
higher than of the typical hard sphere one due to the 
softness of the cages which allows more flexibility in the 
rattling than in the case of rigid hard spheres. 

We now turn to examine the density correlators. In 
Fig. [6] we show the behavior of 4>u(q,t) for large (a) and 
small stars (b) respectively for several wave vectors. As 
for the MSD, the large stars display a marked plateau 
at all wave vectors, lasting for about three decades in 
time. The height of the plateau oscillates as a function of 
the wave- vectors, as commonly found in standard glasses. 
We extract the non-ergodicity parameter of the system 
by fitting large stars 4>n{q,t) with a stretched exponen- 
tial law (Eq. 16 ). The results are shown in Fig.[7]together 
with the corresponding MCT predictions. We find that 
/11(g) oscillates in phase with the static structure fac- 
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FIG. 6: Density correlators for large (a) and small (b) stars 
in the single class for different wave numbers q (see legends). 
The parameters are the same as in Fig. [5] The dashed lines 
are stretched exponential fits (Eq. ( |16[ ) to the curves. The 
inset (b) shows a magnification of the fits for g-values where 
small stars also display a two-step decay. 



tor (not shown), but in contrast to the MCT predictions, 
it shows a marked increase at low q, while beyond the 
first peak the agreement is very good. We find that the 
stretching exponent (3 q as a function of qo\ also oscillates 
in phase with fu(q), and tends to approximately 0.6 for 
large q. Indeed, using this value as b exponent for the 



von Schweidler law (Eq. ( 15 1), we find that all curves are 



well described in the region of departure from the plateau 
(Fig-[6ja)). This result for b agrees well with MCT pre- 
dictions reported above. The estimates for fn{q) from 
the von Schweidler fits are in perfect agreement with that 
obtained from the stretched exponential fits (Fig. [7]) . 

For the small stars correlators (Fig.[6^b)), we find that 
the plateaus are sometimes so small that no accurate 
analysis can be performed. Nonetheless at small wave 
vectors, the stretched exponential fits work quite well 
(see inset). We find finite values /22(<z) just close to the 
first peak of Su(q), suggesting a confinement of the small 
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FIG. 7: Comparison of partial non-ergodicity parameters 
from MCT and simulation in the single glass. The values 
from the simulations are obtained by stretched exponential 



fits (Eq. (16l) for both types of stars. Also shown are results 
from von Schweidler fit (Eq. ( 15 1) for large stars with b = 0.6. 



The simulation parameters are the same as in Fig. [5] The 
two component MCT results refer to the critical parameters: 
0.345, p 2 o\ = 0.12 and 5 = 0.1. 



stars in the frozen matrix of the large stars, which in- 
duces a small loss of ergodicity. Thus a coupling in time 
between large and small stars exists only at the peak of 
Sn(q), which might explain why two-component MCT is 
still capable of correctly describing the general behavior 
of the mixture in this regime. The reason why the modes 
around q « 5af 1 do not couple with modes at different q 
values is due to the fact that the structure factor of the 
small stars is structure-less except for this q value and 
close to the hydrodynamic limit (q — > 0). 

Theoretical descriptions of the dynamical arrest of a 
fluid in a frozen matrix have recently been developed 
[Ml [70] which corresponds to the single glass state found 
here. Fig.|7]shows that /22(g) is in good agreement with 
the MCT predictions, but we do not detect any increase 
of non-ergodicity at large length-scales. While MCT un- 
derestimates the non-ergodicity parameter of the large 
stars at low q, it overestimates the same for the small 
stars compared to the simulation results. 

We can conclude from the combined study of MCT 
and simulations of the single glass that it corresponds 
to a glass formed by the large stars, whose properties 
are those typical of ultra-soft star polymer glasses. The 
presence of the small stars induces a weakening of the sta- 
bility of the glass, until eventually melting it. We observe 
this in the behavior of the large stars' diffusion coefficient 
which grows with increasing density of the additives. The 
small stars are truly ergodic and fully mobile within the 
voids of the glassy matrix. They are never trapped in 
a metastable glass state. Nevertheless, they show a loss 
of ergodicity for length-scales of the order of the static 
structure factor first peak for the large stars. We re- 
mark that this behavior is seen both in binary MCT and 
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FIG. 8: Mean squared displacements for large and small stars 
in the double glass. The densities are p\o\ = 0.41 and pic\ = 
0.1; the size ratio is <5 = 0.4. The horizontal lines indicate the 
estimates of the squared localization length for large and small 
stars, respectively, defined as the inflection point of the MSD. 



simulations, while it could not have been detected by an 
effective one-component treatment of large stars only. 



B. Region II: Double glass 

We now turn to examine the arrested state for low pi 
at large S. We consider the case P\<j\ = 0.41, S = 0.4 and 
pio\ = 0.1, so that the composition of the mixture is the 
same as for the single glass examined above. However, 
the additives are larger and, consequently, also their mass 
is larger. Here the mass ratio is m^/mi = 8.5 • 10~ 2 , 
which is one order of magnitude larger than in the single 
glass. Indeed, the simulations confirm an arrest of both 
species. 

In Fig. [8] we show the MSD for both components in the 
mixture. Differently from the single case, we observe a 
significant slowing down also in the small component be- 
havior, suggested by the intermediate time regime where 
the slope of (r 2 ) versus time (in the double-logarithmic 
plot) is smaller than 1, indicating a subdiffusive regime 
in the dynamics. In the same time window, the large 
stars are more confined than the small ones; their local- 
ization length ^ MD remains equal to that of the single 
glass. For the small stars we can deduce instead a local- 
ization length of the order of 0.5ci, about twice as large 
as that of the large stars, again in agreement with MCT 
predictions. 

Looking at the partial density correlators for the dou- 
ble glass (Fig.[9^a)), we see a close similarity for the large 
stars correlators to those of the single glass. Only the 
time duration of the plateau changes, due to the differ- 
ent relative position of the chosen state point with respect 
to the glass transition. Indeed, it can be seen in Figs, [I] 
and [4] that this state point is 'more liquid' than the one 
we studied for the single glass. Similarly to the previ- 
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FIG. 9: Density correlators for large (a) and small (b) stars 
in the double glass for different values of wave number q (see 
legend). The parameters are the same as in Fig. [8] The 
dashed lines for the large stars are stretched exponential fits 
(Eq. p6|), while for the small stars we use logarithmic fits 
(Eq. |l7P) with t/tmd = 5. 



ous paragraph, we can analyze the decay in terms of a 
stretched exponential and of a von Schweidler law and the 
results are totally equivalent to those of the single glass 
concerning dynamics of the large stars. The resulting 
/11(g) is reported in Fig. 10 together with corresponding 



critical MCT predictions in this regime. 

A completely different picture arises for the small stars 
(Fig. |9jb)). The density correlators cannot be described 
in terms of stretched exponentials. Indeed, they are best 
fitted with second order polynomial in ln(t) of Eq. (17 1 
with t/tmd = 5. In other systems, this anomalous dy- 
namical behavior has been attributed to the competition 
between two glassy states [55l [59j [60l [61], [71]. This seems 
to suggest the possible existence of a higher order MCT 
singularity |63j . although we do not detect an increase 
of the exponent parameter A in the MCT calculations. 
Interestingly, the logarithmic behavior is found only for 
the small stars, while the behavior of the large stars re- 
mains close to a standard liquid-glass transition or A2 
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FIG. 10: Comparison of partial non-ergodicity parameters 
from MCT and simulation in the double glass. The val- 
ues from the simulations are obtained by stretched exponen- 
tial (Eq. (16l) fits for the large stars and from logarithmic 
(Eq. \17\ ) fits for the small stars. The simulation parame- 
ters are the same as in Fig. H] The MCT results refer to the 
critical parameters: p\o\ = 0.345, p2<J\ = 0.05 and S = 0.4. 
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FIG. 11: Mean squared displacements for large and small 
stars in the asymmetric glass. The densities are p\a\ = 0.345 
and p\o\ = 3.4; the size ratio is S — 0.4. he horizontal lines 
indicate the estimates of the squared localization length for 
large and small stars, respectively, defined as the inflection 
point of the MSD. 



singularity. 



Performing the fits following Eq. (17 1, we extract a 



non-ergodicity parameter (Fig. 10 1 for the small compo- 



nent in good agreement with the one calculated within 
MCT close to the liquid-glass transition (P2&1 = 0.05, 
see Fig. [lj. Also, by looking at the (not shown) 
extracted from the fits, we can find the wave vector q* 
at which the decay is purely logarithmic [53J |75] , which 
indicates the dominant length scale in the arrest pro- 
cess, and indeed it was estimated before for the attrac- 
tive glass [55]. Here, h q ~ for qo\ w 4, corresponding 
to a length-scale slightly larger than the nearest-neighbor 
length of the large stars, indicating that small stars are 
mostly trapped in the voids between large stars. Around 
this g*-value there is the expected crossover between con- 
cave and convex shape for ^22(9, t) [S3]- As for the single 
glass case, the estimates of fu(q) from the simulations 
differ from the MCT ones especially at low q, though the 
discrepancy is smaller in this case. 

To summarize, the simulation results are in agreement 
with the MCT predictions discussed above, which suggest 
that the glassy properties of the large stars are identical 
in the single and double glass, both in terms of non- 
ergodicity parameter and localization length, while the 
small stars display a crossover from ergodic to arrested 
anomalous behavior. 



C. Region III: Asymmetric Glass 

We analyze the state point pi = 0.345, p2<jf — 3.4 
and S = 0.4 (5378 particles in the simulation), where the 
diffusivity decreases upon addition of small stars, as can 
be seen in Fig. [4] 



From the MSD reported in Fig. 11 we see that both 
components exhibit a slowing down in the dynamics, 
therefore we have evidence of a second state in which 
both components become arrested, i.e. a second double 
glass. However, the localization lengths that can be ex- 
tracted from the MSD are smaller than those for the dou- 
ble glass discussed above. Indeed, we find Imd ~ 0.14/ci 
for the large stars, corresponding to about half of the cage 
length that we found both in single and double glass, 
while Imd ~ 0.2/ui for the small stars, again at least 
twice as small as that of the other (conventional) double 
glass. These values are again in remarkable agreement 
with those predicted by MCT. The additives do not only 
lead to a reduced localization length, but also deform the 
cages of the big particles, leading to an asymmetric glass 

02- 

The analysis of the density correlators provides evi- 
dence that the dynamics in this double asymmetric glass 
is very different from that of the standard double glass, 
see Fig. 12 Indeed, we cannot fit the correlators of the 
large stars with stretched exponentials, but rather they 
arc found to follow the logarithmic decay of Eq. ( 17 1 with 
t/tmd = 1- On the other hand, the small stars 4>22(l,t) 
are well described by a stretched exponential decay for 
small wave vectors. The logarithmic behavior for large 
stars suggest the possibility of the presence of a higher 
order MCT singularity for the large stars or, at least, a 
competition between two different glassy states. Indeed, 
upon increasing pi at fixed S = 0.4, the large stars are 
passing from the standard glass to a much more local- 
ized one. The same is true for the small stars in this 
region, so that a similar log-behavior is also expected for 
the small stars at some region of the phase diagram, but 
within the present study we are not able to detect that. 
However, we observe a deviation from standard stretched 
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FIG. 12: Density correlators for large (a) and small (b) stars 
in the asymmetric for different wave numbers q (see legend). 
The parameters are the same as in Fig. |11| Here the large 
stars are described by polynomial in ln(t) (Eq. (17l), for the 
small ones we use stretched exponentials (Eq. (16 1). 



exponential behavior at large q also for the small stars. 

From the logarithmic fits for the large stars we de- 
termine the non-ergodicity parameters for large stars, 
again in quite good agreement with MCT predictions (see 
Fig. 131. Also, we can determine the q* where the time 
dependence is logarithmic, i.e. q*a± ~ 16.6, thus invoking 
smaller length scales than the nearest-neighbor distance 
as responsible for the mechanism of arrest. We notice 
that this length is much larger than that found for at- 
tractive glasses in CP mixtures |55j . consistent with the 
fact that here no attractive bonding is present. For the 
small stars, the stretched exponential fits are good at 
small wave-vectors, but become more uncertain for large 
q values. However, the fits show a satisfying agreement 
for the non-ergodicity parameters with the MCT predic- 
tions for the critical parameters p\ — 0.345, pi = 1.27 
and 6 = 0.4. We have also tried to use a Von Schweidler 
exponent b = 0.5, quite close to the MCT predictions, 
which also fits quite well the small stars correlators. 

To recapitulate, also in the asymmetric glass we find 
good agreement between simulation results and MCT 
predictions. Moreover, we find again anomalous dynam- 
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FIG. 13: Comparison of partial non-ergodicity parameters 
from MCT and simulation in the asymmetric glass. The val- 
ues from the simulations are obtained by stretched exponen- 
tial ( Eq. ( 16 1 ) fits for the small stars and from logarithmic 
(Eq. (17l) fits for the large stars with t/tmd = 1- The simu- 
lation parameters are the same as in Fig. [llj The MCT results 
refer to the critical parameters: picrf = 0.345, P2&1 = 1-27 
and S = 0.4. 



ics. However, with respect to the double glass the prop- 
erties of the small and large stars density correlators are 
reversed, i.e. now the large stars can be described in terms 
of stretched exponentials while the small ones exhibit a 
logarithmic decay. This fact, together with the growth 
of A for this case as reported in the MCT results section, 
points to the possible existence of a nearby higher-order 
singularity for the large stars. 



D. Region IV: Attractive Glass and Phase 
Separation 

Finally, we turn to the region of small S and large p2 
where we expect the occurrence of the attractive glass. 
However, this is difficult to observe due to competition 
with phase separation in the mixture. As already dis- 
cussed before, we can only detect an attractive glass with 
binary MCT, because we are not able to observe a bend- 
ing of the ideal glass line in the one-component MCT 
treatment before the integral equation calculations cease 
to converge. 

Nevertheless, in the simulations, we do observe a de- 
crease of the diffusion coefficient before the onset of phase 
separation (Fig. 14 1. After an initial increase of diffu- 
sivity the dynamics becomes slower with increasing pi. 
However, this decrease is not sufficient to observe the 
emergence of a plateau in the MSD and to extract values 
for the nonergodicity parameter from the simulation. 

We then used the static structure factor from the 
(monodisperse) simulations as input for the mode- 
coupling equations and detected a glassy state both in bi- 
nary and one-component MCT for p\cr\ — 0.345, S = 0.1, 
and p2<J i = 20. The partial structure factors are shown in 
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FIG. 14: Diffusion coefficient of the large stars, normalized 
by the value for the reference one-component system Do, for 
p\o\ — 0.345 and 8 — 0.1 as a function of p2- As an inset we 
show the MSD of the large stars for different concentrations 
of small ones. 



Fig. 15 'a). The small wave-length behavior of the struc- 
ture factors shows the typical increase of a state point 
close to phase separation. In Fig. |l5[b) we show the cor- 
responding non-ergodicity parameters whose shape looks 
similar to the one found in previous studies of systems 
with attractive interactions [29 . It is interesting to note 
that the one- and two-component MCT give almost the 
same results for fu(q). Hence, the decrease of the diffu- 
sivity combined with large-scale fluctuations seems to in- 
dicate the possibility of an "arrested phase separation" [5] 
in this case, the interplay beween glass transition and 
phase separation induces the dynamical arrest of the 
denser phase. This mechanism has been recently pro- 
posed to describe dynamical arrest at low and interme- 
diate densities in CP mixtures [73] and other attractive 
systems [74] . 



VII. CONCLUSIONS 

In conclusion, we have identified different glassy states 
in a binary mixture of star polymers by using extensive 
computer simulations and mode coupling theory. In the 
plane spanned by size asymmetry between the two com- 
ponents and concentration of the added small star, there 
is a liquid lake which reveals multiple reentrant behav- 
iors. Both the upper and lower threshold concentrations 
for the liquid-glass transition behave nonmonotonically 
with the size asymmetry 5. At high concentrations pi 
of small stars and S = 0.4, we have identified a new, 
asymmetric double glass type which is characterized by 
a strong localization of the small particles and anisotropic 
cages. This is different both from the usual double glass 
at low concentrations pi and from the hard sphere dou- 
ble glass [33] . Recent experiments confirm our theoretical 
findings [42] , At low S and high pi there is a competi- 




FIG. 15: (a) Partial structure factors p\o\ = 0.345, 5 = 0.1, 
and p2C"i = 20. Note the small-q behavior of all three curves, 
which signals phase separation, (b) The corresponding non- 
ergodicity parameters predicted by MCT. 



tion between vitrification and phase separation. In all 
the explored phase space, we find remarkable agreement 
between simulation results and MCT predictions. 

An interesting open question is which minimum arm 
numbers /i and /2 are necessary for the formation of a 
high-p2 glass, as we know that one-component star poly- 
mer systems with less than ss 40 — 45 arms do not ever 
glassify |34j . An additional interesting question concerns 
the presence of higher order singularities in the wide pa- 
rameter space that can be explored with star-polymer bi- 
nary mixtures. Indeed, the MD results provide evidence 
of a non-conventional slowing down of the dynamics in 
the proximity of the double and of the asymmetric glass. 
It is important to find out if MCT predicts any sharp 
glass-glass transition between single and double glass, or 
between double and asymmetric glass. This will be the 
subject of a future study. A more exhaustive theoreti- 
cal study is requested to properly map the MCT results 
with the numerical ones and provide a guide-line for lo- 
cating, if existing, the higher order singularities, similarly 
to what has been done for the case of short-range attrac- 
tive colloids [551 172"]. 
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